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ABSTRACT 

We use ID radiative transfer simulations to study the evolution of H2 gas-phase (H~ 
catalyzed) formation and photo-dissociation regions in the primordial universe. We find 
a new positive feedback mechanism capable of producing shells of H2 in the intergalactic 
medium, which are optically thick in some Lyman- Werner bands. While these shells 
exist, this feedback effect is important in reducing the H2 dissociating background flux 
and the size of photo-dissociation spheres around each luminous object. The maximum 
background opacity of the IGM in the H2 Lyman- Werner bands is th 2 ~ 1 — 2 for a relic 
molecular fraction xh 2 = 2 x 1CP 6 , about 6 times greater than found by Haiman et al. 
(2000). Therefore, the relic molecular hydrogen can decrease the photo-dissociation rate 
by about an order of magnitude. The problem is relevant to the formation of small pri- 
mordial galaxies with masses Mum i$ 10 8 M , that rely on molecular hydrogen cooling 
to collapse. Alternatively, the universe may have remained dark for several hundred 
million years after the birth of the first stars, until galaxies with virial temperature 
T v ir > 10 4 K formed. 

Subject headings: Galaxies: formation - Cosmology: theory 



1. Introduction 

Many recent observations point to a flat geometry and a significant cold dark matter (CDM) 
content. Recent measurements of the Cosmic Microwave Background (CMB) anisotropies (de 
Bernardis et al. 2000; Lange et al. 2001; Tegmark et al. 2000), large-scale structure (Hamilton & 
Tegmark 2000), Lya power spectrum (Croft et al. 1999; McDonald et al. 2000), and distances of 
high redshift SNe (Perlmutter et al. 1998; Riess et al. 1998; Garnavich et al. 1998) are all consistent 
with a flat CDM model of the universe, in which about 1/3 the total energy density of the universe 
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is dark matter (DM) and 2/3 is "vacuum energy". Big-Bang nucleosynthesis fits to the observed 
D/H abundances suggest a small baryon density, with closure parameter Q^h 2 ~ 0.02 (Buries &; 
Tytler 1998). Thus, the baryons are a small fraction of the DM. In the linear phase of structure 
formation, the baryons simply trace the DM perturbations on large scales. 

In CDM cosmologies, small mass DM virialized objects form first. Larger halos form later 
from the merger of smaller subunits (hierarchical scenario). Since the gas temperature prior to 
reionization is about 100 K, the gas can collapse in objects with Mdm > 10 4 — 10 5 M . Small mass 
DM halos will then produce the first luminous objects (Pop III) if the baryons collapsed into the 
DM potential wells and heated to the virial temperature can dissipate their internal energy support. 
The lack of metals in the primordial universe makes molecular hydrogen the only coolant below 
T ~ 10 4 K. Therefore, unless some H2 is formed, the baryons can not cool in shallow DM potential 
wells with Mdm ^ 10 s [(1 + z)/10]~ 3 / 2 M , corresponding to a virial temperature T„ r < 10 4 K. 
The formation of H2 occurs mainly through the chemical reaction, 

H + H-^H 2 + e~. (1) 

Together with a similar reaction (i?^ - + H — > H2 + H + ), this is the only possibility for virialized 
objects with T v i r < 10 4 K to cool and eventually form stars (Lcpp & Shull 1984). Despite the fact 
that the physics is much simpler, theoretical models of the formation and fate of Pop III objects 
must confront the lack of observational constraints and uncertainties on initial conditions. The fate 
of the Pop III objects is still unresolved or at least controversial. The main complications of the 
model are the feedback processes between luminous objects and the formation and destruction of 
H2 that can affect star formation on both local (interstellar medium, ISM) and large (intergalactic 
medium, IGM) scales. 

Molecular hydrogen is photo-dissociated by photons with energies between 11.1 and 13.6 eV 
(Lyman- Werner bands) through the two step Solomon process (Stecher & Williams 1967). Thus, 
the radiation from the first stars (or quasars) could dissociate H2 before the IGM is reionized 
{negative feedback). The H formation process instead relies on the abundance of free electrons. A 
positive feedback is therefore possible in a gas partially ionized by an X-ray background or direct 
flux of Lyman continuum radiation. Depending on the relative importance of the positive/negative 
feedback, star formation could be triggered or suppressed after the formation of the first stars. 

Understanding the star formation history at high redshift is important for modeling the metal 
enrichment of the IGM. It may also influence IGM reionization if, as pointed out in Ricotti & Shull 
(2000), the contribution to the Lyman continuum (Lyc) emissivity in the IGM at high redshift 
is dominated by low-mass galactic objects. The contribution to the emissivity from these objects 
is important because their spatial abundance is higher and because the fraction, (f eS c), of Lyc 
radiation that escapes from a low-mass galaxy halo is larger than for more massive objects. 

Tegmark et al. (1997) used a simple collapse criterion, that the gas cooling time must be less 
than the Hubble time, to determine the minimum H2 fraction a virialized object must have in 
order to collapse. This minimum mass is a decreasing function of the virialization redshift and it is 
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somewhat dependent on the cooling function used. Tegmark et al. (1997) found a minimum mass of 
Mom ~ 10 6 M at z ~ 30 using the Hollenbach & McKee (1979) molecular cooling function. Abel 
et al. (1998) and Fuller k, Couchman (2000) instead found a minimum mass an order of magnitude 
smaller using Lepp & Shull (1984) and Galli & Palla (1998) cooling functions respectively. Objects 
of this masses can form around redshift z ~ 30 from 3<r perturbations of the initial Gaussian density 
field. More realistic models (Omukai & Nishi 1998; Nakamura & Umemura 1999) confirm this basic 
result. 

Abel et al. (2000) and Bromm et al. (1999) used high-resolution numerical simulations to 
study the formation of the first stars. Abel et al. (2000) used 3D adaptive mesh refinement (AMR) 
simulation with dynamic range about 2 x 10 5 to follow the collapse of the first objects starting 
from cosmological initial conditions (standard CDM cosmology). They found an approximatively 
spherical protogalaxy of 10 6 M with a collapsing core of ~ 100 M . Similar results have been 
found by Bromm et al. (1999) using SPH simulations starting from ad hoc initial conditions. For 
some initial conditions and larger halo masses they could form disk protogalaxies with more than 
one star formation region. Unfortunately, without taking into account radiative transfer effects and 
SNe explosions, it was not possible to determine the mass or the initial mass function (IMF) of the 
first stars. Local effects of UV radiation on star formation have been studied with semi-analytical 
models (Omukai &: Nishi 1999; Glover & Brand 2001). The main result is that massive stars are 
effective in suppressing cooling in small protogalaxies, but star formation can continue in larger 
galaxies. If the ISM is clumpy, star formation is not suppressed in the denser clumps. 

The photo-dissociating UV radiation can affect the H2 abundance over large distances if the 
IGM is optically thin in the Lyman- Werner bands. Haiman Sz Loeb (1997) and Haiman et al. 
(2000) computed self-consistently the rise in the dissociating UV background (UVB) using the 
Press-Schechter formalism for sources of radiation. They found that H2 in the IGM had a negligible 
effect on the build-up of the UVB because of its small optical depth th 2 ^0.1 and because the 
photo-dissociation regions around the first luminous objects overlap at an early stage. Therefore, 
the negative feedback of the UVB suppresses star formation in small objects. However, if the first 
objects are mini-quasars, the produced X-ray background is strong enough to cancel out the UVB 
negative feedback, and reionization could be produced by small objects. Contrary to Haiman et al. 
(2000) result, Ciardi et al. (2000) found no negative feedback from the UVB. They estimated the 
mean specific intensity of the continuum background at z = 20 in the 912 — 1120 A Lyman- Werner 
bands, J u ~ 10~ 27 erg s _1 cm -2 Hz" 1 sr _1 . This value is four orders of magnitude smaller than 
the value found by Haiman et al. (2000). Even if the photo-dissociation regions overlap by z = 25, 
only at z < 15 does the UVB become the major source of radiative feedback. At higher redshifts, 
direct flux from neighboring objects dominates the local photodissociation rate (Ciardi et al. 2000). 
Machacek et al. (2001) included the feedback effect of the photo-dissociating UV background in 3D 
AMR simulations of the formation and collapse of primordial protogalaxies. They used a box of 1 
Mpc 3 comoving volume, and resolution in the maximum refinement zones of about 1 pc comoving. 
Their results confirm that a photo-dissociating background flux of Flw ^5x 10 -24 erg s _1 cm -2 
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Hz 1 delays the cooling and star formation in objects with masses in the range 10 5 - 10 7 M . They 
also provided an analytical expression for the collapse mass threshold, for a range of UV fluxes. 

The simple observation that star formation in our Galaxy can trigger further star formation 
in a chain-like process (McCray & Kafatos 1987) suggests that positive feedback effects could be 
important also at high redshift. Haiman et al. (1996) found that UV irradiation of a halo can 
enhance the formation of H2 and favor the collapse of the protogalaxies. Ferrara (1998) estimated 
that the formation of H2 behind shocks produced during the blow-away process of Pop III can 
replenish the relic H2 destroyed inside the photo-dissociation regions around those objects. Finally, 
it has been noticed that positive feedback is possible if the first objects emit enough X-rays (Haiman 
et al. 2000; Oh 2001; Venkatesan et al. 2001). 

In this paper, we explore a new positive feedback effect that could be important in regulating 
star formation in the first galaxies. Ahead of the ionization front, a thick shell of several kpc of 
molecular hydrogen (with abundance xh 2 = n H 2 / n H 10 ) forms because of the enhanced elec- 
tron fraction abundance in the transition region between the H II region and the neutral IGM. This 
shell, which we call a positive feedback region (PFR), can be optically thick in the H2 Lyman- Werner 
bands, depending on the source turn-on redshift Zi, luminosity, and escape fraction (f esc ) of ioniz- 
ing radiation. We find that the photo-dissociation spheres around single objects are smaller than 
in previous calculations where the effects of the PFRs and H2 self-shielding have been neglected. 
PFRs could also be important in calculating the build-up of the dissociating background. A self- 
consistent calculation of feedback effects on star formation, including the new feedback process 
discussed in this study, will be the subject of a subsequent paper. We will use a 3D cosmological 
simulation with radiative transfer, in order to understand what regulates the star formation at high 
redshift. 

This paper is organized as follows. In § 2 we describe chemical and cooling/heating processes 
included in the radiative transfer code. In § 3 we discuss some representative simulations and find 
analytical formulae that fit the simulation results. In § 4 we use this analytical relationship to 
constraint the parameter space where negative or positive feedback occurs. In § 5 we discuss the 
opacity of the IGM to the photo-dissociating background. We summarize our results in § 6. 

Throughout this paper, when we talk about "small halos", we refer to objects with virial 
temperatures T v i r < 10 4 K that rely on the presence of molecular hydrogen in order to collapse. 
Conversely, "large halos" are objects with T v i r > 10 4 K that can cool down and collapse even if 
H2 is completely depleted. The cosmological model we adopt is a flat ACDM model with density 
parameters (f&A,f2o>^&) = (0.7,0.3,0.04), and Hubble constant Hq = 70 km s -1 Mpc -1 . Thus, the 
baryon density is 0&/i 2 = 0.0196, compatible with D/H inferences from Big-Bang nucleosynthesis 
(Buries & Tytler 1998). 
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2. 



Radiative Transfer in the Primordial Universe 



We use ID radiative transfer simulations to investigate the evolution of H2 photo-dissociation 
regions in the primordial universe. We are interested in studying the evolution of the abun- 
dances and ionization states of hydrogen, helium, and molecular hydrogen around a single ioniz- 
ing/dissociating source that turns on at redshift z%. Hubble expansion and Compton heating/cooling 
are also included. 



The radiative transfer equation in an expanding universe in comoving coordinates has the 
following form: 



where J V {R) [erg s _1 cm~ 2 Hz -1 sr _1 ] is the specific intensity and a = (1 + z em )/(l + z a b s ), with 
z em and z a bs the emission and absorption redshifts, and e u and k v the emissivity and absorption 
coefficients. If the scale of the problem is smaller than the horizon size, c/H, and \dJ v /d\rvv\ ~ J u 
(true in the case of continuum radiation), we can neglect the cosmological terms. In § 2.4 we 
show that we can also use the classical radiative transfer equation for lines if we account for the 
cosmological redshift of the radiation (or line cross sections) . Finally, if the emission and absorption 
coefficients do not change on time scales shorter than a light-crossing time, (L/c), we can solve at 
each time step the stationary classical radiative transfer equation. For the case of a single point 
source, the solution is the attenuation equation 



2.1. ID Non-equilibrium Radiative Transfer 



l dJ v {R) l dJ„(R) H ( dJ u (R) 
c dt a dR c V d\nv 





MR) = 



Mo) 



exp[-T v (R)], 



(3) 



R 2 



where 




(4) 



with rii the number density and <jj the absorption cross section of species i. 



2.2. Chemistry and Initial Conditions 



We solve the time dependent equations for the photo-chemical formation/destruction of 8 
chemical species (H, H + , H _ , H2, ELj~, He, He + , He ++ ), including the 37 main processes relevant to 
determine their abundances (Shapiro h Kang 1987). We use ionization cross sections from Hui & 
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Gnedin (1997) and photo-dissociation cross sections from Abel et al. (1997). The energy equation 
is 



where e g = (3kT/2)riH (1 + x(He) + x e ), with nu the hydrogen number density, and x(He) and x e 
the helium and electron fractions respectively. Note that nu is a function of time because of the 
Hubble expansion. Also, x e is time dependent. The cooling function A includes H and He line and 
continuum cooling (Shapiro & Kang 1987) and molecular ro-vibrational cooling excited by collisions 
with H and H2 (Martin et al. 1996; Galli & Palla 1998). As mentioned above, adiabatic cosmic 
expansion cooling is also included. The heating term T includes Compton heating/cooling and 
photoionization/dissociation heating. We solve the system of ODEs for the abundances and energy 
equations, using a 4 th -order Runge-Kutta solver. We switch to a semi-implicit solver (Gnedin & 
Gnedin 1998) when it is more efficient (i.e., when the abundances in the grid are close to their 
equilibrium values). Convergence analysis showed that a logarithmic grid with 400 cells in space 
and 400 cells in frequency (for the continuum flux) coordinates is sufficient to achieve convergence 
within a 10% error. The spectral range of the radiation is between 0.7 eV and 1 keV. 

The primordial helium mass fraction is Yp = p(He)/pb = 0.24, where pb is the baryon density, 
so that x(He) = Yp/4(1 — Yp) = 0.0789. The initial values at z = Z{ for the temperature and 
species abundances in the IGM are: T = 10 K, xh 2 = 2 x 10~ 6 , and x e ~ x H + = W~ 5 / (hQ^n^ 2 ) = 
6.73 x 10 -4 . The initial abundance of the other ions is set to zero. We explore two cases: (i) sources 
embedded in the constant density IGM and, (ii) sources inside a virialized halo. For case (ii) we 
use initial conditions and a density profile to match the numerical simulation of Abel et al. (2000). 

The density distribution of massive stars in the halo is crucial for determining {f esc ), the escape 
fraction of ionizing photons (Ricotti & Shull 2000). If all the stars are located at the center of the 
DM potential, it is trivial to find (f eS c) for a given baryonic density profile, because the calculation 
is reduced to the case of a single Stromgren sphere, 



where a H = 2.59 x 10 -13 cm 3 s _1 is the case-B recombination rate coefficient at T = 10 4 K and 
Jin (cm -3 ) is the hydrogen number density. Assuming that the gas profile is given by solving the 
hydrostatic equilibrium equation in a DM density profile (Navarro et al. 1997), and that the Lyc 
photon luminosity, So [photon s _1 ], is proportional to the baryonic mass of the galaxy, we find, 




(6) 




,3 



if(l + z)<1.22(e//,) 1 /3, 
if (l + z)>1.22(e// 9 ) 1 /3. 



(7) 



Here, z is the collapse redshift of the halo, S = (1.14 x 10 49 s~ l )ef g (M DM /10 6 M ), e is the star 
formation efficiency normalized to the Milky Way, and f g is the collapsed gas fraction; see Ricotti 
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Sz Shull (2000) for details. We expect f g to be small since the DM potential wells of this objects are 
too shallow to hold photo-ionized gas. When (/ esc ) is not zero, the number of photons absorbed 
is such that all the gas in the halo is kept ionized. Since we assume spherical symmetry, we can 
only place the source at the center of the halo. From the steep rise of (/ eS c) m equation (7), it 
is clear that (f esc ) is essentially either 1 or for the majority of the redshifts and star formation 
efficiencies. In a more realistic geometry, (f esc ) turns out to be small but not zero. We can make 
a more realistic calculation if we place the grid outside the virial radius in the constant-density 
IGM and reduce the ionizing flux by a factor (fesc) 1 - The escape fraction from the halo can be 
estimated with more accurate calculations (Dove et al. 2000; Ricotti & Shull 2000) or left as a free 
parameter. 

2.3. Photon Production 

The spectral energy distribution (SED) of the source is assumed to be either a constant power 
law, F v oc v~ a , with a = 1.8 for the case of a mini-quasar (QSO), a constant (at t = 0) SED from 
a population of metal- free stars for the Pop III source (Tumlinson & Shull 2000), or an evolving 
SED of stars with metallicity Z = 0.001 (Leitherer et al. 1995) for the Pop II source. We consider 
the two extreme cases of an instantaneous burst of star formation or continuous star formation. 

In Figure 1 we show the SED of a zero metallicity population for instantaneous star formation, 
in which 100 M of gas is converted into stars with a Salpeter IMF, a lower mass cut-off at M\ ow = 1 
Mq and upper cut-off at M up = 100 M . The solid line shows the non-evolved SED (at t = 0). 
We also show the photo- ionization and photo-dissociation cross sections for H (dotted line), H2 
(dashed line), and (dot-dashed line). The H2 photo-dissociation cross section in the Lyman- 
Werner bands is the average cross section of the lines between 11.2 and 13.6 eV according to Abel 
et al. (1997). In Figure 2 we show the SED of the Z = 0.001 metallicity population with the same 
parameters as for the Pop III in Figure 1. 

We have explored the effects of changing the IMF and mass cut-off. The main difference is 
that the Lyc luminosity So can change by as much as an order of magnitude for a steeper IMF 
(a = 3.3) or for a Salpeter IMF (a = 2.35) with M up = 30 Mq. The Lyc luminosity at t = 
generated by converting instantaneously 10 6 M of baryons into stars is So ~ 10 53 [photon s _1 ] 
for the Pop III SED and S ~ 6 x 10 52 [photon s" 1 ] for the Pop II SED. The specific flux from the 
object is F v = L u /AttR 2 [erg s _1 cm~ 2 Hz -1 ] where L v is the specific luminosity. In a previous 
conference proceedings (Ricotti et al. 2001) we erroneously underestimated the specific intensity 
by a factor 4-7T for the Pop II and Pop III SEDs. Therefore the stated photon luminosities in those 
figures should be a factor of 4ir bigger. 



1 For the sake of simplicity we simply reduce the ionizing flux by a factor (f eS c) without changing the SED of the 
source that, in a realistic case, should become harder. Indeed, as noted in § 3, in the simulations where the source is 
at the center of a halo the emerging ionizing photons are both reduced and harder 
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1 10 1 10 2 

hi/[eV] 

Fig. 1. — SED of a zero metallicity population for instantaneous star formation, in which 100 Mq 
of gas is converted into stars with a Salpeter IMF, a lower mass cut-off at M[ ow = 1 M Q and upper 
cut-off at M up = 100 M . The solid line shows the non-evolved SED (at t = 0). We also show the 
photo- ionization and photo-dissociation cross sections for H (dotted line), H2 (dashed line), and 
(dot-dashed line). The H2 photo-dissociation cross section in the Lyman- Werner bands is the 
average cross section of the lines between 11.2 and 13.6 eV according to Abel et al. (1997). 
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Fig. 2. — SED of the Z = 0.001 metallicity population with the same parameters as for the Pop III 
in Figure 1. The lines show the SED at t = 1, 5, 10, 100 Myr after the burst (Leitherer et al. 1995). 
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2.4. Radiative Transfer for the Lyman- Werner Bands 

In Figure 3 we show a representative simulation after t ~ 16 Myr from a burst of a Pop III 
objects turning on at Z{ = 19 with total Lye photon luminosity So = 1.2 x 10 49 s _1 . We find a shell 
of H2 formation just in front of the H II region that we call a positive feedback region (PFR), where, 
in some cases, the abundance of H2 is much higher that the relic abundance of xh 2 = 2 x 10 -6 . 
Figure 4 shows the gas temperature as a function of the comoving distance from the source at, 
t ~ 0.05, 1, and 22 Myr for the same object shown in Figure 3. Even if the temperature is T ~ 10 4 
K in the PFR, collisional dissociation of H2 is not very effective at the typical densities of the IGM 
or in the outer part of galaxy halos. 

In other cases, we find that the column density of H2 formed ahead of the ionization region is 
high enough (N > 10 14 cm -2 ) that some Lyman- Werner bands are optically thick at line center. 
In order to properly include this effect, we solve the radiative transfer, not only for the continuum 
radiation but also in lines. Since the H2 and H I lines are very narrow (because the IGM temperature 
is T ~ 10 K [Auth/v ~ 10~ 6 ]), we need extremely high spectral resolution in the frequency range 
of 10.2 — 13.6 eV. Fortunately, we only need to resolve 76 absorption lines longward of 912 A [R(0), 
R(l), P(l) Lyman lines and R(0), Q(l), R(l) Werner lines] arising from the J=0 and J=l rotational 
levels of v = in the H2 ground electronic state. Thus, we only need high spectral resolution in the 
vicinity of the absorption lines. This limits the number of the frequency bins to just 3 x 10 4 instead 
of ~ 10 6 for the uniform sampling. At densities uh % 10 cm~ 3 , the population of upper rotational 
or vibrational states of H2 is negligible even at T = 10 4 K. We also include some H I Lyman-series 
resonance lines that are close in frequency to H2 lines. For these lines we use a Voigt profile, since 
they are usually in the damping-wing regime of the curve of growth. For the Lyman- Werner lines 
we use Gaussian profiles since the optical depth does not exceed a few. We redshift the specific 
flux according to the distance from the source, so that the maximum redshift of a line depends 
on the grid physical size. For our simulation the maximum frequency redshift is smaller than the 
distance of two consecutive Lyman or Werner line "triplets". Therefore we need to consider only 
a limited frequency range around each triplet. We resolve the frequency profile of each one of 
the 76 absorption lines in the Lyman- Werner bands with about 100 grid points if T ~ 10 4 K and 
with a few points if T ~ 10 K. Another complication arises from the fact that our spatial grid 
is logarithmic. Therefore, at a large enough distance from the source, the line can be redshifted 
several line widths inside each grid cell. We will show in § 5 that, for the simple case of constant 
density, temperature, and molecular abundance, the optical depth can be calculated analytically. 
The line profile, in this case, is similar to a top-hat function with r = XH 2 n H(z)o'f osc ^Xi/H(z) if 
Vi — Av < v < Vi and zero otherwise, where Av = H(z)Ax/\ (see Figure 11). Here, H(z) is 
the Hubble constant, a = (ire 2 /m e c) is the classical cross section, f osc ,i is the oscillator strength, 
and Xi is the wavelength of the i th line. In Figure 5 we show an example of the Lyman- Werner 
bands spectrum emerging from the PFR after t = 50 Myr from the source turn-on. The source is a 
Pop III object with Lyc photon luminosity So ~ 10 49 s" 1 turning on at z = 30. Each panel shows 
the spectrum at progressively higher resolution. 
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R [kpc] 

Fig. 3. — Abundances as a function of the comoving distance from the source at t ~ 100 Myr for a 
Pop III object turning on at Zj = 30 with Lyc photon luminosity Sq = 1.2 x 10 49 s _1 . The H and 
He I ionization fronts are at R ~ 5 kpc and just ahead of the fronts is the PDF region with an H2 
abundance xh 2 — 5 x 10 -4 . The thick (thin) line shows xh 2 including (excluding) line radiative 
transfer in H2 Lyman- Werner bands. 
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Fig. 4. — Gas temperature as a function of the comoving distance from the source at times t ~ 
0.05, 1, and 22 Myr after a Pop III object turns on at Zi = 19 with Lyc photon luminosity 
S = 1.2 x 10 49 s" 1 . 
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Fig. 5. — Emerging spectrum in the Lyman- Werner bands at t = 50 Myr after a Pop III object 
turns on at Zi = 30 with Lyc photon luminosity Sq ~ 10 49 s _1 . Each panel shows the spectrum at 
progressively higher resolution. We only show the Lyman series absorption lines of hydrogen that 
are important in reducing the photo-dissociation rate of H2. The molecular absorption lines are 
produced by the optically thick PFR. The solid line shows the region of the spectrum in which we 
calculate line radiative transfer. The dotted line shows the SED of the source. 
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3. Results 

3.1. Simulation Results 

Preceding the ionization front, a thick shell of several kpc of molecular hydrogen (with xh 2 
up to 10~ 4 ) forms because of the enhanced electron fraction in the transition region from the 
H II region to the neutral IGM. This shell can be optically thick in some Lyman- Werner H2 lines 
depending on the redshift, source luminosity, and escape fraction. This has two main consequences: 
(1) the photo-dissociation fronts around the sources slow down and their final radii are smaller than 
in the optically thin case; (2) optically thick PFRs could reduce the intensity of the cosmological 
background in the Lyman- Werner bands. 

In Figures 6-8 we show the isocontours of \og(xH 2 ) as a function of time and comoving distance 
from the source. The thick lines show the analytical fits for the position of the photo-dissociation, 
formation and ionization fronts as a function of time. Figure 6 (top) shows an example where 
positive feedback is unimportant. The emitting object has a Pop II SED, So ~ 10 50 s _1 , and a 
turn-on redshift Z{ = 30. The source is at the center of a baryonic halo with density profile n oc R~ 2 . 
Figure 6 (bottom), instead, shows an example where positive feedback is important. The source 
has Pop III SED, So ~ 10 50 s -1 , and a turn-on redshift = 20. The photo-dissociation front moves 
so slowly that the PFR shell gets close to it, producing an enhancement of H2 abundance instead of 
a depletion. In Figure 7 (top) we show a QSO SED, So — 10 51 s" 1 , and a turn-on redshift Zi = 20. 
The mini-quasar produces fronts similar to the Pop III object. In Figure 7 (bottom) the source 
is a Pop II object with instantaneous burst of star formation. After t ~ 10 7 yr, when OB stars 
start to explode as supernovae, the H II region recombines, triggering the formation of H2. Finally, 
Figure 8 (bottom) shows the effect of including radiative transfer in the Lyman- Werner lines; the 
photo-dissociation front slows down considerably with respect to the optically thin case (top). 

The UVB background affects these results in a simple way. The H2 abundance is xh 2 oc 
1/(F S + Flw) where F s and F^w are the fluxes in the Lyman- Werner band from the source and 
from the background respectively. In the optically thin case F s ~ 8.4 x lO^-R^So/lO 49 s- 1 ) 
[erg s _1 cm~ 2 Hz -1 ]. Therefore, the background flux equals the source flux at a distance from the 
source R crit ~ (8.4 x 10" 22 /^LH/) 1/2 (5'o/10 49 s" 1 ) 1 / 2 kpc. At a distance R = 3R crit the background 
flux will reduce our calculations of xh 2 by a factor of 10. 

3.2. Analytical Fits 

Let us introduce some relevant time scales. The collisional recombination time scales for neutral 
hydrogen and protons are tni = t rec (l — x e )/x 2 and tun = t rec /x e respectively, where x e is the 
electron fraction number density. Here we define t rec = 1/ {a^ nn) where = 2.59 x 10 -13 cm 3 
s _1 is the case-B recombination coefficient at T = 10 K and nn(zi) = (1.7 x 10~ 7 cm~ 3 )(l + Zj) 3 is 
the hydrogen number density for O5/1 2 = 0.019. The age of the universe tn and t rec at the redshift 
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Fig. 6. — Contours of log(xjj 2 ) as a function of time and comoving distance from the source. The 
thick lines show the analytical fits of the photo-dissociation front (D), ionization front (I) and 
formation region (F). (Top) example of negative feedback: Pop II SED, Sq ~ 10 50 s _1 , turn-on 
redshift Zi = 29. The source is at the center of baryonic halo with density profile n oc R~ 2 . 
(Bottom) example of positive feedback: Pop III SED, Sq ~ 10 50 s _1 , turn-on redshift Zi = 19. 
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Fig. 7.— Same as Figure 6. (Top) QSO SED, So ~ 10 51 s _1 , turn-on redshift Zi = 19. (Bottom) 
fossil H II Region: Pop II SED (instantaneous star formation law), Sq ~ 10 49 s -1 , turn-on redshift 
Zi = 14. 
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Fig. 8. — Same as Figure 6. (Top) Pop III SED, So ~ 10 49 s -1 , turn-on redshift Zi = 30. (Bottom) 
same but including line radiative transfer in H2 Lyman- Werner bands. 
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Zi when the source turns on are: 



-3/2 



1 7 20 (±±3.) Myr. (9) 



If the density of the IGM is constant and uniform, the radius Rj of the ionization front is, 



Rs 



1 — exp 



t 

ire( 



l 

(io) 



where Rs = (3Sot rec / 'Aimnizi)] 1 ^ 3 is the Stromgren radius and t is the time the source is on. The 
typical time for the H II region to reach the Stromgren radius in a non-expanding IGM equals 
the H + recombination time. A source turning on at redshifts Z{ < 18, for which A = tnjt r( , c < 2, 
produces an H II region that never reaches the Stromgren radius (Donahue & Shull 1987; Shapiro 
& Giroux 1987). In comoving coordinates, it expands forever because of cosmic dilution. This is 
true only if the source continues to shine with constant photon luminosity Sq. We therefore remind 
the reader of another typical time scale, namely the time t^urst — 5 — 20 Myr during which an 
instantaneous burst of star formation produces a substantial amount of ionizing photons. 

Based on the simulations, we find analytical relationships that are good fits to the comoving 
front radii (radius of the ionization sphere, Ri, radius of the photo-dissociation sphere, Rp>, and 
outer radius of the PFR, Rf), the peak H2 abundance xh 2 ,f, comoving thickness, ARp and column 
density, Np, of the PFR as a function of time. We also estimate the time scale, tp, on which a 
region in the IGM lies inside the PFR. These relations, without taking into account line radiative 
transfer, are: 

R D = (i6kpc)(^-)[pS 0A9 f 2 (t)]h, (11) 



10 



Rj = (21kpc)[5 0) 4s</ cac )/i(t)]3 ) (12) 
ARj = (5 x 10- 5 7 kpc) (1 + f' 3 (S 0A9 (fesc)t) - 18 (13) 

= (3.2 7 Mpc) (1±^)°' 7 (g °' 49 g^ t)a18 , (14) 

AR F = (15) 

where So, 49 is the ionizing photon luminosity in units of 10 49 [photon s -1 ], (3 and 7 are fitting 
parameters (Table 1) that depend on the SED of the object (/? = 7 = 1 for a Pop III SED), and t 
is the time the source was on in Myr. ARj is the comoving thickness of the ionization front. The 
functions f\ and fi are given by 

hit) = i8 '^*lp X exp [£^{t c E 2 {\/t c ) - E 2 {\)) Myr, (16) 

/ 2 (t) = 3t H [l-t-V 3 ] Myr, (17) 
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with t c = 1 + t/tfj, A = tn/trec an d E n (x) is the exponential integral of order n (Donahue & 
Shull 1987; Shapiro & Giroux 1987). Note that ~ t if t < t H and f 2 (t) ~ t if A < 2 (i.e. 

2j < 18). Table 1 shows the normalization parameters f3, 7 and the ionizing photon luminosity 
<So,49 produced by an instantaneous burst of ~ 250 M© of gas into stars for the SEDs from Pop III, 
Pop II (continuous star formation) and mini-quasars. 

Here the photo-dissociation front is defined as the locus where the molecular abundance is 
xjj 2 = 10~ 6 , half the relic H2 fraction. In the optically thin regime, that is a good approximation 
in this case, the profile of the H2 abundance inside the dissociation sphere is exponential: xh 2 = 
2 x 10 -6 exp[— ln(2)i?£)/r]. The photo-dissociation front slows down at t > tn, after the initial 
i 1 / 2 expansion law. At late times the dissociation front approaches the maximum radius Rd = 
(340a) kpc[(l + ^)/10] 1 / 4 (/3S" ,49) 1/2 , where a = [1 - (1 + Zi)' 1 / 2 ] 1 / 2 . The relation (1 + z) = 
(1 + Zi)t c 2 ^ 3 is used to relate time and redshift. The maximum thickness of the ionization front 
ARj = (7 7 Mpc)[(l + Zi)/10] a43 (S* ,49 (fesc)) - 18 is reached at t = 0M H - 

When the ionization front is inside the halo with a density profile n oc R~ 2 , where R is the 
radius, we find log.Rj oc (logt) 2 . The emerging spectrum in this case is harder. The photo- 
dissociation front, if we do not take into account self-shielding in H2 lines, is independent of the 
chosen density profile of the gas. The width of the PFR inside the halo is smaller by a factor re -1 , 
but its column density increases as n. 

We have made some numerical "experiments" to test which species have abundances equal to 
their equilibrium values. Only H~ can be safely considered in equilibrium. Nevertheless, for H2 the 
equilibrium abundance is not a good approximation, it is an helpful estimate since can be derived 
analytically (Donahue & Shull 1991; Abel et al. 1997) and will give us the functional form for the 
fitting formula. If we consider only the formation process through H~, with rate kf{T)xtf- and 
photo-dissociation by photons with hu > 0.755 eV, the equilibrium H2 abundance in the PFR is 
xh 2 ,f of n H{z)kf{T)x H -F~ l oc (l+z)Rj /(fiSo^g), where F s oc /3So,4g(l+z) 2 Rj 2 is the dissociating 
flux. In all the simulations, after the PFR shell starts to build up, the abundance of H~ reaches an 
equilibrium value x H - ~ 5 x 10 -7 . The best fit differs slightly from the expression above because 



Table 1. Fit parameters for the analytical formulae in § 3.2 



Object 


-50,49 


P 


7 




(250 M burst) 






Pop II 


2.8 


2.4 


0.67 


Pop III 


3.2 


1.0 


1.0 


Quasar 




0.8 


1.17 
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of a redshift dependence of the formation rates: 

x H2 , F = (I.74xl0~ 9 ) (1 + Z)L3 ^ (18) 

P>J0,49 

= i.exio-^i^y^-^-lffc) 2 /^. (19) 

The maximum value of xh 2 ,f, reached at t = (10/3)i#, is given by 

xffi? = 5.9 x 10- 4 (^y 3 S~^(f esc )l. (20) 

From equations (11-19) it is easy to derive useful relationships to quantify the importance of 
positive feedback as a function of the unknown free parameters. We will use these relationships in 
the next section. The ratio 

= 0.76 (I±*) {f esc )- l/ ^ l/2 sl%g{t) (21) 

reaches the maximum value of 

/ d \ max / 1 i \ 3/4 

(f^) =l-6(^) {fes C )~^^Sl% (22) 

at t ~ 0M H where the function g(t) = [/i(i)] 1/2 [/2(i)] _1/3 scales as - * 1/6 if * < *if and A < L 
The function has a maximum value of g = 0.75i]^ 6 = 2.1[(1 + z^/lO]" 1 / 4 at i = 0.9tf/. 

It is also useful to estimate the H2 column density, Nf, of the PFR, and the time scale, £f, on 
which a region in the IGM can be engulfed by the PFR: 

UXh^ARf , 1Rvln i 3 - 2 ,7 + ((/esc)*)°- 85 ^ 



Nf ax = (9 x 10 14 cm" 2 )] f l±^l\ LS (/egc) °' 85 ( 24 ) 



with the maximum value of 

.mi x m 14 ,.,„ -, , _ 

at i ~ 0.9iff, and 

Finally, we derive the ratio R D /R F = Rd/{Ri + ARf) = {Rd/Ri){1 + ARf/Ri)' 1 that gives 
us an idea of the fraction of the IGM volume in which the molecular hydrogen is destroyed: 

-1 

(26) 



Rd_Rd 
R f ~ Ri 



I. + 51' i^-) ' 7 (SoMfesc)t)-°- 15 



10 
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with the maximum value of 



Rp 
R F 



) 



max 



0.15/3 



a/2 r 




n -0.18 



7 



(fesc) ( 



,1/3 
'0,49 



(27) 



at t ~ 0.9<h- 



4. Discussion: Negative or Positive Feedback ? 



In this section, we use the analytical relationships found in § 3.2 to quantify the importance 
of the positive feedback as a function of the free parameters of the model. The free parameters 
are: the Lyc escape fraction, {f e sc}, the star formation efficiency e normalized to that in the Milky 
Way, and the collapsed gas fraction, f g . For a fixed cosmology, the ionizing photon luminosity is 
So = (1-14 x 10 49 s- 1 )e/ 3 (M z?M /10 6 M ) [see Ricotti & Shull (2000) for details]. We quantify the 
importance of the positive feedback by showing in Figure 9 isocontours of constant ratio, Rd/Rf, 
of the photo-dissociation front radius to the formation front radius 10 Myr after the source turned 
on. Here the radius of the dissociation front Rd is defined as the locus where the relic molecular 
abundance has dropped to xh 2 = 10 -8 , therefore it is smaller by a factor ln(200)/ln(2) ~ 7.6 with 
respect to the values given in the previous paragraph. Clearly, if Rd/Rf < 1 the dissociation 
region does not exist. Instead, the H2 abundance could increase with respect to the relic value 
xh 2 ~ 2 x 10~ 6 (as an example see Figure 6 [bottom]). Therefore, Rd/Rf gives us an idea of 
the mean molecular abundance in the IGM and the opacity of the IGM to the photo-dissociating 



We find that the existence of positive feedback depends crucially on {f esc ) and the SED of the 
first objects (Pop III, Pop II or AGN). Figure 9 shows the parameter space where positive feedback 
is possible in the particular case of a Pop III SED. The contour lines show Rd/Rf as a function of 
redshift and mass of the halo. The shaded region shows the parameter space where Rd/Rf < 1, in 
which case the PFR fills up the region between the H II region and the photo-dissociation front. At 
small redshifts, the boundary of the shaded region is produced by the additional constraint that the 
H2 abundance of the PFR has to be xjj 2 > 10~ 5 . The choice of this value is somewhat arbitrary, 
but the main purpose is to show the H2 abundance in the PDF as a function of the redshift and 
halo mass. The thick solid line at the bottom is the minimum mass to collapse as a function of 
redshift according to the criteria of Abel et al. (1998). The two dashed lines show objects with 
virial temperatures of T v i r = 10 3 K and T v i r = 10 4 K; protogalaxies that lie above T v i r = 10 4 K 
are not subject to negative/positive feedback because they can cool by H I line excitation (Lya). 
The minimum mass for collapse is calculated assuming that the object lies outside of a photo- 
dissociation region or a PFR and that the dissociating background is negligible. The amount of H2 
formed in a just-virialized object does not depend on the initial H2 abundance, but it is sensitive 
to the electron fraction. For this reason, protogalaxies less massive than those found in the Abel 
et al. (1998) calculations can form inside a PFR. Finally, we note that the presence of PFRs will 
increase the opacity of the IGM in the Lyman- Werner lines, further decreasing the intensity of the 



background. 
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Redshift (z) 




Redshift (z) 



Fig. 9. — The contours show Rd/Rf for our standard model and cosmological parameters (§ 1). 
The shaded region (Rd/Rf < 1) shows the redshift and halo masses where positive feedback is 
dominant. The thick solid line shows the minimum mass needed to collapse according to Abel 
et al. (1998). Parameters: (f eS c) = escape fraction, e = star formation efficiency normalized to the 
Milky-Way, f g = collapsed gas fraction, So = (1.14 x 10 49 s~ 1 )e/ 9 (M DM /10 6 M ). The dot-dashed 
lines show the collapse redshift of la and 3a perturbations according to linear theory. (Top) Pop III 
SED, e=l,f g = l and (f esc ) = 0.2. (Bottom) Pop III SED, e = 1, f g = 0.1 and (f esc ) given by 
the analytical formula derived in Ricotti & Shull (2000). 
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dissociating UV background. All these effects have to be included in a 3D cosmological simulation 
with radiative transfer in order to understand their global effect on structure formation and the 
star formation history. 



5. IGM optical depth in the Lyman- Werner bands 

In this section we estimate the IGM optical depth in the Lyman- Werner bands that arises from 
the absorption in 76 narrow H2 lines from J = 0,lint> = 0at different redshift. This calculation is 
analogous to the one in Haiman et al. (2000). We repeat it because we disagree with their technical 
analysis of H2 line scattering and fluorescence. Haiman et al. (2000) assumed that only the fraction 
of absorptions that decay to the dissociating continuum remove photons from the UVB (about 
11% of the time). For the other 89% of the absorptions, they assume that the Lyman- Werner 
photon is re-emitted at the same frequency, with no net effect. This is not correct. Following 
the absorption in the Lyman or Werner bands, the molecule decays to a variety of ro- vibrational 
levels in the ground electronic state because the quadrupole transitions within the electronic excited 
state are much slower. The H2 either decays to the 6 3 S+ (anti-bonding) state which decays to the 
vibrational continuum (pre-dissociation) or to a bound vibrational-rotational level in the electronic 
ground state (X 1 ^^). At moderate UV intensities, the subsequent infrared cascade through the 
bound levels of the ground electronic state is entirely determined by the radiative decay rates (Black 
k, Dalgarno 1976). From the cascade probabilities, Shull (1978) computed that a Lyman- Werner 
photon has a 14% probability, on average, to be re-emitted at the same frequency. Approximately 
12% of the absorption transitions dissociate and 74% fluoresce to excited (v, J) levels of X l Y,+ . 
Roughly speaking, the probability to decay to one of the other 14 bound vibrational levels is about 
6%. Therefore, only 14% of the time the Lyman- Werner photon is resonantly scattered following an 
absorption. The other 86% of the time, H2 dissociates or the photon is split into infrared and less 
energetic (about 1 eV energy loss) UV photons that are removed from the Lyman- Werner bands 
right away or after a few absorptions. 

Therefore, the H2 optical depth of the IGM is given by 

7re 2 n ^> rz em dt 

Tu(z b,Zem) = 22 foscji 1 ~ fi,v"=o) dz' C— n H (z')x H2 {z')(t){v' , Vi) (28) 

m eC . J z ob az 

where z Q b and z em are the observer and emission redshifts, v' = u{l + z em )/(l + z b), tih{z) = 
(1.12 x 10~ 5 cm" 3 )fi fe /i 2 (l - Y P )(1 + z) 3 is the neutral H number density in the IGM, <j> is the line 
profile, fosc,i is the oscillator strength, and fiy=o the probability to decay to the ground vibrational 
level of the X 1 ^^ state for the i th line calculated from Black Sz Dalgarno (1976). The maximum 
redshift interval a UV photon can travel before it is absorbed by a neutral atom corresponds to the 
redshift between two H I resonances in the higher Lyman series, the "dark screen" approximation 
of Haiman et al. (2000). Therefore a photon is subject to the absorption from a subset of the 76 
Lyman- Werner lines, n\(v) < i < ri2(^), where n\(i>) is the first Lyman- Werner line with frequency 
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just above v, and n-iiy) > n\(v) is the Lyman- Werner line with frequency just below the next 
higher H I Lyman line with frequency > v. 

It is possible to write an approximate analytical solution of equation (28), if we assume Gaus- 



3 x 10- 7 ViT}£ M is the 



sian line profiles, 4>g(v, v i) = (I/v^^j) ex P[ — ( u ~ v i) 2 1 ^IXi where Az^ 
Doppler width of the H2 line i, or Lorentzian line profiles 4>l{ v i v i) = (^i/^ 7r )/[( u ~ v i) 2 + (Fi/2) 2 ], 
where Tj = 7i/27r and 7* = ^2iA(i — > I) is the natural width of the i th H I line. For a constant 
xh 2 = 2 x 10~ 6 molecular abundance, we have 



/i+ \ 3 / 2 ™ 2< ^ 
r^zot, * em ) ^ 3.28 x 10 15 ( ) ^ / MC)i (l - f iy , =0 )v] /2 <S>(p' 



(29) 



where, 



$(1/) = < 



1 

V 7T 



-f(^f)-erf(^) 

(^f) - arctan (^) 



arctan 



for Gaussian profile 
, for Lorenzian profile 



(30) 



and erf(x) is the error function. In Figure 10, we show the H2 total optical depth (i.e., up to the 
maximum z em visible to the observer) of the IGM in the Lyman- Werner bands at redshifts z = 30 
(solid line) and z = 15 (dashed line). The opacity at energies less than Ly/3 {hv < 12.09 eV) 
is produced by the H2 Lyman lines. At energies higher than Ly/3, the H2 Werner lines are also 
important. The maximum opacity that we find is r ~ 2, about 6 times higher than that found by 
Haiman et al. (2000). The background flux is thus reduced by an order of magnitude if we assume 
an average molecular fraction xh 2 = 2 x 10~ 6 . 

For the Gaussian line profile, the expression inside the square brackets of equation (30) tends 
to the asymptotic value of 2 when (1 + z em )/(l + z Q b) > A^. The maximum optical depth of the 
i th line is therefore 



,{z ob ) ~ XH 2 n H (z = 0)—af o 
no 



(1 + z ob )^ 



~ 0.1 



2 x 10" 



12 eV 

hvj 



foSCy 

1(F 



1 + z ob 
30 



3/2 



(31) 



In Figure 11, we show the optical depth through a constant density and molecular fraction gas for 
a single line. When the source redshift z em increases, the profile becomes wider and the optical 
depth tends to the asymptotic value r max in equation (31). 

Finally, we calculate the total optical depth of the IGM as a function of distance from a 
single source. We also include in the calculation the Lyman series lines of H I and we investigate 
the importance of the damping-wings of the heavily saturated lines. The mean optical depth for 
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hu [eV] 



Fig. 10. — Total H2 optical depth of the IGM in the Lyman- Werner bands. The solid line is for an 
observer at redshift z a b = 30 and the dashed line at z a b = 15. We assume a constant xh 2 = 2 x 10 -6 . 
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Fig. 11. — In the plot we show the function [erf(x') — erf(x)], where x' = [y' — v^/Avi and 
x = [y — fi)/Avi. According to equation (29), this function is proportional to the optical depth of 
an H2 line through a gas of constant density and molecular fraction. When the source redshift z em 
increases with respect to the observer redshift z Q b, the profile becomes wider, and the optical depth 
tends toward the asymptotic value T max in equation (31). The profiles are shown at log[(z em — 
Zob)/{l + z ob)] = —3.5, —4, —4.5, —7. The gas temperature is either T = 10 4 K (solid lines) and 
T = 10 K (dashed lines). 
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R [Mpc] (comoving) 

Fig. 12. — Total optical depth, {tigm), for photo-dissociation of H2 in the IGM as a function of the 
comoving distance from a source at z em = 30. {tigm) includes the opacity of H2 and H I lines and 
is weighted on the photo-dissociation rate; see equation (33). The solid line is calculated assuming 
a Voigt profile for the Lyman series hydrogen lines, therefore including the effect of the damping 
wings. The dashed line is calculated neglecting H I damping wings. The insert is a zoom of the 
inner few Mpc, and the arrow indicates the distance at which the absorption lines are redshifted 
about one line width for a gas at T = 10 4 K. 



DRAFT: February 1, 2008 



28 



photo-dissociation, (tjgm), can be calculated by solving the following equations: 

kdiss = const ^ fosc,ifdiss / dv'-^j 4>{y' ,v { ), (32) 
i=i Jv hv 

76 

^ ss exp(-(T JGM » = const ^/ OSCil / dlss / du' ^{v' ^e'^M^ ^ m ,z oh ) (33) 

»=i ^ 

where /c^ ss is the photo-dissociation rate and where const=47ra;// 2 ?7.//(2; f ) )o' with o = (ire 2 /m e c). 
Note that tigm differs from t v in equation (29) since it includes the opacity of both H2 and H I lines. 
In Figure 12 we show (tjgm) as a function of the comoving distance (and z {, < z em ) from a source 
at z em = 30. The solid line is calculated assuming a Voigt profile for the Lyman series hydrogen 
lines, therefore including the effect of line saturation. The dashed line is calculated neglecting line 
saturation. The insert is a zoom of the inner few Mpc, and the arrow indicates the distance at 
which the absorption lines are redshifted about one line width for gas at T = 10 4 K. Both Figure 12 
and Figure 5 show clearly that the opacity at comoving distances from the source of less than 1 
Mpc is produced in the PFR shell. The IGM opacity starts to be important around 10 Mpc. This 
reduces the photo-dissociation radius of the strongest sources at high redshift. 

We remind the reader that all the calculations in this section assume a constant molecular 
fraction in the IGM of xjj 2 = 2 x 10~ 6 . In order to calculate a realistic average molecular fraction 
in the IGM as a function of the redshift, we need to calculate the filling factor of the photo- 
dissociation regions, PFR and the H2 formed inside relic H II regions. Only with a 3D radiative 
transfer simulation will it be possible to treat self-consistently the aforementioned feedback effects. 



6. Discussion and Summary 

The results of this work reinforce the possibility that the population of small mass primordial 
galaxies (Pop III) could exist, without being immediately suppressed by radiative feedback. A 
quantitative answer to this question has to be given by high-resolution cosmological simulations. 
Until the build-up of the dissociating background, local feedback effects regulate star formation. 
When the star formation is suppressed in a halo, it is reasonable to expect that H II regions start to 
recombine. The relic H II region, as shown in § 3, produces new molecular hydrogen and perhaps 
a second burst of star formation. This is especially effective at high redshifts where the density of 
the IGM and protogalaxies is higher and the physical scales are smaller. Thus, at these redshifts, 
it is easier to believe that star formation is bursting rather than continuous. The picture that 
could emerge from numerical simulations is that H II regions and photo-dissociation regions in 
the IGM are short-lived instead of continuously expanding, as in the reionization simulations of 
Gnedin (2000). This will have important effects on the distribution of the metals and the chemical 
evolution of galaxies. 

In summary, our key results are: 
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(1) We have found a new positive feedback effect on the formation of H2. Each source of radiation 
produces a shell (PFR: positive feedback region) of H2 in front of the H II region with a thickness 
of several kpc and peak abundance xh 2 ~ 10 -4 . The H2 column density of the PFR is typically 
Np(H2) ~ 10 14 — 10 15 cm -2 . Fossil H II regions, if they exist, are an important mechanism of H2 
production, both in the IGM and inside protogalaxies. 

(2) The PFR can be optically thick in the lines of the H2 Lyman- Werner bands. The implication 
is twofold: (i) the photo-dissociation region around each single source is about 1.5 times smaller 
than in the optically thin case; (ii) the H2 optical depth of the IGM increases. In a later paper, 
we will calculate the IGM optical depth and the effect on the intensity of the photo-dissociating 
background by means of cosmological simulations. 

(3) We provide analytical formulae that fit the simulation results, in order to make a parametric 
study of the importance of positive feedback as a function of redshift. The most important parame- 
ters are (fesc) an d the SED of the sources (Pop III, Pop II or AGN). If (f esc ) is not extremely small, 
PFRs have important effects on the opacity of the IGM to the H2 photo-dissociating background 
and the size of photo-dissociation regions. 

(4) The background opacity of the IGM in the H2 Lyman- Werner bands is about unity if xh 2 = 
2 x 10~ 6 . Therefore, if the relic molecular hydrogen is not immediately destroyed, it can decrease 
the photo-dissociating background flux by about an order of magnitude. 

The aforementioned results are the foundations on which we will construct a cosmological 
simulation in which these effects are treated self-consistently. 

This work was supported by the Theoretical Astrophysics program at the University of Col- 
orado (NSF grant AST 96-17073 and NASA grant NAG5-7262). We thank Mark Giroux for a 
critical review of the manuscript and Jason Tumlinson for giving us his Pop III SED. 
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